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Abstract. Multi-planet systems detected until now are in most cases characterized by hot- 
Jupiters close to their central star as well as high eccentricities. As a consequence, from a 
dynamical point of view, compact multi-planetary systems form a variety of the general A''- 
body problem (with N ^ whose solutions are not necessarily known. Extrasolar planets are 
up to now found in prograde (i.e. direct) orbital motions about their host star and often in 
mean-motion resonances (MMR). In the present paper, we investigate a theoretical alternative 
suitable for the stability of compact multi-planetary systems. When the outer planet moves on 
a retrograde orbit in MMR with respect to the inner planet, we find that the so-called retro- 
grade resonances present fine and characteristic structures particularly relevant for dynamical 
stability. We show that retrograde resonances and their resources open a family of stabiliz- 
ing mechanisms involving specific behaviors of apsidal precessions. We also point up that for 
particular orbital data, retrograde MMRs may provide more robust stability compared to the 
corresponding prograde MMRs. 
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1. Introduction 

To identify the dynamical state of multi-planetary systems, we use the MEGNO tech- 
nique (the acronym of Mean Exponential Growth factor of Nearby Orbits; Cincotta 
& Simo 2000). This method provides relevant information about the global dynamics 
and the fine structure of the phase space, and yields simultaneously a good estimate of 
the Lyapunov Characteristic Numbers with a comparatively small computational effort. 
From the MEGNO technique, we have built the MIPS package (acronym of Megno Indi- 
cator for Planetary Systems) specially devoted to the study of planetary systems in their 
multi-dimensional space as well as their conditions of dynamical stability. 

Particular planetary systems presented in this paper are only used as initial condition 
sources for theoretical studies of 3-body problems. By convention, the reference system 
is given by the orbital plane of the inner planet at t = 0. Thus, we suppose the orbital 
inclinations and the longitudes of node of the inner (noted 1) and the outer (noted 2) 
planets (which are non-determined parameters from observations) as follows : ii — Q° and 
Vli =0° in such a way that the relative inclination and the relative longitude of nodes 
are defined at f = as follows : v = «2 — *i = ^2 and Q.^ = VI2 — = ^2- The MIPS 
maps presented in this paper have been confirmed by a second global analysis technique 
(Marzari et al. 2006) based on the Frequency Map Analysis (FMA; Laskar 1993). 



2 



Julie Gayon & Eric Bois 



2. Fine structure of retrograde resonance 

Studying conditions of dynamical stability in the neighborhood of the HD 73526 two- 
planet system (period ratio: 2/1, see initial conditions in Table 1), we only find one stable 
and robust island (noted (2)) for a relative inclination of about 180° (see Fig. [5^). Such 
a relative inclination (where in fact ii = 0° and 12 — 180°) may be considered to a 
coplanar system where the planet 2 has a retrograde motion with respect to the planet 1. 
From a kinematic point of view, it amounts to consider a scale change of 180° in relative 
inclinations. Taking into account initial conditions inside the island (2) of Fig. la, we 
show that the presence of a strong mean-motion resonance (MMR) induces clear stability 
zones with a nice V-shape structure, as shown in Fig. lb plotted in the [ai, ei] parameter 
space. Let us note the narrowness of this V-shape, namely only about 0.006 AU wide 
for the inner planet (it is 5 times larger in the Jupiter-Saturn case). A similar V-shape 
structure is obtained in [02, 62] with about 0.015 AU wide. Due to the retrograde motion 
of the outer planet 2, this MMR is a 2:1 retrograde resonance, also noted 2:-l MMR. 
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Figure 1. Panel (a); Stability map in the [ir,r2r] non-determined parameter space of the 
HD 73526 planetary system (see Table 1). Panel(b): Stability map in the [ai, ei] parameter space 
for initial conditions taken in the stable zone (2) of panel (a). Note that masses remain untouched 
whatever the mutual inclinations may be; they are equal to their minimal observational values. 
Black and dark-blue colors indicate stable orbits (< Y >= 2 ± 3% and < Y >= 2 ± 5% respec- 
tively with < Y >, the MEGNO indicator value) while warm colors indicate highly unstable 
orbits. 
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Figure 2. Stability maps in the [ir,r2r] parameter space. Panel (a): initial HD 82943 planetary 
system (see Table 1). Panel (b): scale reduction of the HD 82943 planetary system according to 
a factor 7.5 on semi- major axes. Masses in Panel (a) and Panel (b) are identical. Color scale is 
the same as in Fig. 1. 
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3. Efficiency of retrograde resonances 

Fig. 2 exhibits stability maps in the [ir, ^r] parameter space considering a scale re- 
duction of the HD 82943 planetary system (see Table 1) according to a factor 7.5 on 
semi-major axes (masses remaining untouched). The dynamical behavior of the reduced 
system (Fig. 2b) with respect to the initial one (Fig. 2a) points up the clear robustness 
of retrograde configurations contrary to prograde ones. The "prograde" stable islands 
completely disappear while only the "retrograde" stable island resists, persists and even 
extends more or less. Even for very small semi-major axes and large planetary masses, 
which should a priori easily make a system unstable or chaotic, stability is possible with 
counter-revolving orbits. 

In the case of the 2:1 retrograde resonance, although close approaches happen more 
often (3 for the 2:-l MMR) compared to the 2:1 prograde resonance, the 2:-l MMR 
remains very efficient for stability because of faster close approaches between the planets. 
A more detailed numerical study of retrograde resonances can be found in Gayon & Bois 
(2008). 



Elements 


HD 73526 


HD 82943 


HD 128311 


HD 160691 


HD 202206 


Mstar (Me) 


1.08 ±0.05 


1.15 


0.84 


1.08 ± 0.05 


1.15 


m sin ii (Mj) 


2.9 ±0.2 


1.85 


1.56 ±0.16 


1.67 ±0.11 


17.4 


2.5 ±0.3 


1.84 


3.08 ±0.11 


3.10 ±0.71 


2.44 


a (AU) 


0.66 ± 0.01 


0.75 


1.109 ±0.008 


1.50 ±0.02 


0.83 


1.05 ±0.02 


1.18 


1.735 ± 0.014 


4.17 ±0.07 


2.55 




0.19 ±0.05 


0.38 ±0.01 


0.38 ± 0.08 


0.20 ±0.03 


0.435 ± 0.001 


e 


0.14 ±0.09 


0.18 ±0.04 


0.21 ±0.21 


0.57 ±0.1 


0.267 ±0.021 


uj (deg) 


203 ±9 


124.0 ± 3 


80.1 ± 16 


294 ±9 


161.18 ±0.30 


13 ±76 


237.0 ± 13 


21.6 ±61 


161 ±8 


78.99 ± 6.65 


M (deg) 


86 ±13 





257.6 ±2.7 





105.05 ± 0.48 


82 ±27 


75.21 ± 1.96 


166 ±2 


12.6 ± 11.2 


311.6 ±9.5 



Table 1. Orbital parameters of the HD 73526, HD 82943, HD 128311, HD 160691 and HD 202206 
planetary systems. Data sources come from Tinney et al. (2006), Mayor et al. (2004), Vogt et 
al. (2005), McCarthy et al. (2004) and Correia et al. (2005) respectively. For each system and 
each orbital element, the first line corresponds to the inner planet and the second one to the 
outer planet. 



Data sources 


Period ratio 


Prograde MMR 


Retrograde MMR 


HD 73526 


2/1 


17 


500 


HD 82943 


2/1 


755 


1000 


HD 128311 


2/1 


249 


137 


HD 160691 


5/1 


e 


320 


HD 202206 


5/1 


e 


631 



Table 2. Statistical results. For each type of MMR (prograde or retrograde), 1000 random 
systems have been integrated in the error bars of each data source. The proportion of stable 
systems over 1000 is indicated in each case, s designates a very small value that depends on 
the size of the random system size. Data sources come from Tinney et al. (2006), Mayor et al. 
(2004), Vogt et al. (2005), McCarthy et al. (2004) and Correia et al. (2005) respectively (see 
Table 1). 
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4. Occurrence of stable counter-revolving configurations 

The occurence of stable two-planet systems including counter-revolving orbits appears 
in the neighborhood of a few systems observed in 2:1 or 5:1 MMR. New observations 
frequently induce new determinations of orbital elements. It is the case for the HD 160691 
planetary system given with 2 planets in McCarthy et al. (2004) then with 4 planets in 
Pepe et al. (2007). Hence, systems related to initial conditions used here (see Table 1) have 
to be considered as academic systems. Statistical results for stability of these academic 
systems are presented in Table 2, both in the prograde case (v = 0°) and in the retrograde 
case (v — 180°). For each data source, 1000 random systems taken inside observational 
error bars have been integrated. Among these random systems, the proportion of stable 
systems either with prograde orbits or with counter-revolving orbits is given in Table 2. In 
all cases, a significant number of stable systems is found in retrograde MMR. Moreover, 
in most data sources, retrograde possibilities predominate. 



5. Resources of retrograde resonances 

The 2:1 (prograde) MMRs preserved by synchronous precessions of the apsidal lines 
(ASPs) are from now on well understood (see for instance Lee & Peale 2002, Bois et al. 
2003, Ji et al. 2003, Ferraz-Mello et al. 2005). The MMR- ASP combination is often very 
effective; however, ASPs may also exist alone for stability of planetary systems. Related 
to subtle relations between the eccentricity of the inner orbit (ei ) and the relative apsidal 
longitude A(D (i.e. (Di —uJ2): Fig. 3 permits to observe how the 2:1 retrograde MMR brings 
out its resources in the [Aw, ei] parameter space : 

• In the island (1) (i.e. inside the [a, e] V-shape of Fig. lb), the 2:-l MMR is combined 
with a uniformly prograde ASP (both planets precess on average at the same rate and 
in the same prograde direction). 

• In the island (2) (i.e. outside but close to the [a, e] V-shape of Fig. lb), the 2:-l 
Tiear-MMR is combined with a particular apsidal behavior that we have called a rocking 
ASP (see Gayon & Bois 2008): both planets precess at the same rate but in opposite 
directions. 

• The [Aw, ei] map also exposes a third island (3) that proves to be a wholly chaotic 
zone on long term integrations. 

Let us note that the division between islands (1) and (2) is related to the degree of 
closeness to the 2:-l MMR. 




Figure 3. Stability map in the [A(D, ei] pa- 
rameter space. A similar distribution of stable 
islands is obtained in [Aw, 62]. Color scale and 
initial conditions are the same as in Fig. 1 with 
in addition the ir and Q.^ values chosen in the 
island (2) of Fig. la. 
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6. Conclusion 

We have found that retrograde resonances present fine and characteristic structures 
particularly relevant for dynamical stability. We have also shown that in cases of very 
compact systems obtained by scale reduction, only the "retrograde" stable islands sur- 
vive. From our statistical approach and the scale reduction experiment, we have expressed 
the efficiency for stability of retrograde resonances. Such an efficiency can be understood 
by very fast close approaches between the planets although they are in greater number. 

We plan to present an Hamiltonian approach of retrograde MMRs in a forthcoming 
paper (Gayon, Bois, & SchoU, 2008). Besides, in Gayon & Bois (2008), we propose two 
mechanisms of formation for systems harboring counter-revolving orbits. Free-floating 
planets or the Slingshot model might indeed explain the origin of such planetary systems. 

In the end, we may conclude that retrograde resonances prove to be a feasible stabilizing 
mechanism. 
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